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We demonstrate fractal noise in the quantum evolution of wave packets moving either ballistically 
or diffusively in periodic and quasiperiodic tight-binding lattices, respectively. For the ballistic case 
with various initial superpositions we obtain a space-time self-affine fractal ty(x,t) which verify the 
predictions by Berry for "a particle in a box" , in addition to quantum revivals. For the diffusive case 
self-similar fractal evolution is also obtained. These universal fractal features of quantum theory 
might be useful in the field of quantum information, for creating efficient quantum algorithms, and 
can possibly be detectable in scattering from nanostructures. 



I. INTRODUCTION 



A quantum wave packet ^(x,t) evolving from a uni- 
form initial superposition state $?(x,0) displays interest- 
ing quantum interference phenomena which involve self- 
affine fractal quantum evolution of the probability den- 
sities and quantum revivals. This was shown by Berry 
[1] in a study of quantum waves confined in boxes by 
Dirichlet boundary conditions. It was pointed out that 
the choice of hard wall boundaries introduces quantum 
fractal noise in the time evolution of uniform determin- 
istic systems obeying the Schrodinger's equation, which 
is expressed via a self-affine fractal probability density 
|\&(a;, t)| 2 . Moreover, over long periods of time at integer 
multiples of the revival time T r , the wave packet ^(x,t) 
returns and reconstructs in its initial form [2] . In a wide 
class of circumstances the quantum revivals are almost 
perfect with both space and time periodic \^(x, t)\ 2 . 

Most studies are concerned with spectra having a 
quadratic dependence on the quantum number, such as 
the simple case of the infinite potential well Ej = Eij 2 , 
j = 1,2,..., where perfect and fractional revivals occur 
at integer and fractional multiples of the revival time T r , 
respectively [1]. For quadratic spectra with wave packets 
localized in energy space around a quantum number j 
the classical time T c i = nh/(Exj) and the much longer 
revival time T r = wh/Ei are the only time scales [3]. 
For more complex energy dispersions, such as the Ry- 
dberg states Ej = Ei/j 2 j = 1,2,... higher order time 
scales exist apart from the classical, e.g. 1st order, 2nd 
order, etc, recurrences, as shown for the "return to the 
origin" probability in the book by Peres [4]. The clas- 
sical revolution time for a wave packet starting around 
the jth orbit is T c i — irfrj /(—E\), while the 1st order 
recurrences occur at integer multiples of the revival time 
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T r = nhj /(—3Ei) [2-4]. The classical time scale T c i de- 
fines the time up to which the system follows a classical 
orbit while the revival time scale T r resembles the Ehren- 
fest time where the classical Liouville density ceases to 
describe the system, since quantum interference makes 
itself obvious with gradual spreading and distortion of 



the wave packet. The classical and revival times are in- 
versely proportional to the energy eigenvalue separation 
and its variation, respectively [4] . The full quantum limit 
t — » oo is approached for much longer times. 

The quantum evolution of wave packets which con- 
sist of a large number of simultaneously excited eigen- 
states is related via the uncertainty principle to ultra- 
short laser pulses of femtosecond duration. In this way 
one can study the evolution of a coherent superposition 
of many quantum states. The quantum evolution of such 
superpositions can be realised in the physics of Rydberg 
atoms [5], molecules [6], semiconductor quantum wells 
[7], mode propagation of Bose-Einstein condensates in 
magnetic wave guides [8], laser prepared states in ionic 
traps [9], etc. On the other hand, from studying squeezed 
wave packets with reduced spatial uncertainty, in order to 
behave like particles as much as possible, one gains a deep 
insight into the semi-classical limit where quantum and 
classical physics meet according to the correspondence 
principle. The quantum evolution is qualitatively milder 
than the classical one and it remembers the initial state, 
due to the reversible and unitary nature of the quantum 
evolution transformation exp(—iHt/h). For classical evo- 
lution the initial state always converges to a stationary 
distribution. We show that quantum memory effects ex- 
ist and a self-affine fractal space and time distribution of 
the probability density \^(x, t)\ 2 . Quantum revivals also 
appear in the "return to the origin" probability. They 
might be interesting for fundamental issues, such as the 
use of quantum mechanics in computation [10]. It was 
shown that space-time interference which repeats itself 
at discrete time intervals can create a subatomic quan- 
tum counter [11] and can help towards implementation of 
factorization for an integer number into primes exploring 
fractional revivals [12] (see also the JV-slit interferometer 
in [13]). 

Our study may also be relevant to another aspect of 
quantum information theory, which is the newborn sub- 
ject of quantum random walks (e.g. the Hadamard walk) 
[14]. A quantum walk is iterated by combining the usual 
translation with a rotation in "coin" -space, where a quan- 
tum coin-flip replaces the classical coin-flip. Of course 
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one should avoid any measurement at intermediate times 
[15]. The quantum walk has an improved behavior over 
the classical random walk since it is characterized by bal- 
listic (quadratic dependence of the mean square displace- 
ment in time), instead of the usual diffusive (mean square 
displacement linear in time) of the classical evolution, ft 
is a general feature of one-particle quantum algorithms 
to speed up quadratically when compared to classical dy- 
namics [14]. 

The quantum evolution is analogous to the classical 
phenomenon of geometrical optics known as the Talbot 
effect [2]. Talbot inspected a coarsely ruled diffraction 
grating illuminated with white light through a magnify- 
ing lens. When the lens held close to the grating, the 
rulings appeared in sharp focus, as expected. However, 
moving the lens away, the rulings of the grating first 
appeared blurred, as expected, but then reappeared in 
sharp focus at multiples of a particular distance, known 
as the Talbot distance. In the Talbot effect the pattern of 
intensity is doubly periodic, both across the grating and 
as a function of distance from the grating. More discus- 
sion on the Talbot effect and its analogies with the quan- 
tum problem can be found in [2] . Although the quantum 
revivals are semiclassical in nature quantum interference 
exists since these semiclassical phenomena appear only 
close to the quantum limit t — > oo. 

In this paper, as an alternative to continuous media 
we aim to study quantum evolution of wave packets mov- 
ing in tight-binding lattices. A discrete and finite space 
is important for simulations on a finite computer with 
the possibility of mapping it to a quantum computer. 
This study is a general scheme for discetizing space when 
one deals with linear equations (such as the quantum 
Schrodinger's equation), which turns out to be particu- 
larly useful when the geometry plays an important role. 
It complies with the necessity to exploit possibilities for 
the information space, which is also discrete, opposite 
to the continuous classical phase space. This convenient 
computational tool for implementing space is favored by 
the presence of an underlying lattice in solid state appli- 
cations. We discuss two cases: (i) the ballistic motion 
in the case of zero site potential where quantum inter- 
ference from wave scattering at the hard walls creates a 
sort of quantum "randomness" similarly to that of the 
"particle in a box" , and (ii) a diffusive case via a special 
quasi-periodic potential of critical strength [16]. It must 
be stressed that quantum diffusive evolution is very dif- 
ferent from classical diffusion although the space fractal 
dimension of the quantum probability is the same. The 
spectrum is of the cos type and the stationary eigenstates 
are simple sin elementary trigonometric functions so that 
the evolution can be characterized by an infinity of other 
time scales apart from T c i, T r . In the ballistic case the 
motion remains ballistic for all t since the tight binding 
lattice has no small t semiclassical limit. 



Our results confirm the presence of fractal fluctuations 
of the probability distribution during the ballistic and 
diffusive quantum evolution, in agreement with Berry 
[1]. Moreover, for all kinds of quantum noise the frac- 
tal noise is described by the predicted universal fractal 
dimensions D x , D t and D x=t of the self-affme curves. 
For the ballistic motion the mean square distance in- 
creases proportionally to the square of the elapsed time 
while for quantum diffusion the mean square distance in- 
creases proportionally to time. The self-affinity seen in 
the noisy quantum evolution is familiar in non-linear sys- 
tems. Therefore, techniques from non-linear physics can 
be used to understand some of the mysteries of quan- 
tum mechanics, particularly quantum evolution in open 
systems in the presence of time-dependent potentials. 



II. BALLISTIC QUANTUM EVOLUTION 
A. Space-time structures 



We have created space-time structures on the discrete 
orthonormal finite lattice space basis {\x),x = 1,2, ...,N} 
for various continuous times t. For the ballistic quan- 
tum evolution in a line the mean-square displacement 
grows quadratically (x 2 ) = 2t 2 . It can be expressed 
via the stationary cigcnsolutions of the Hamiltonian H 
through the evolution operator exp(—iHt/fi), since H 
is time-independent. The static Schrodinger's equation 
H\j) = Ej\j), j = 1, 2, N corresponding to the matrix 
H , which is tridiagonal for a one-dimensional lattice with 
zero diagonal entries and constant off-diagonal elements, 
has stationary eigenvalues, eigenvectors Ej = 2 cos -^q-j- , 

b) = \J ivTT 12x=i sin jtjk \ x }' The quantum evolution 
of the initial ket |^(0)) can be expressed as 

|*(t)) =e- im / n \V(0)) 

N 

= £ e " iB ' t/R IM'l*(o)> (i) 

i=i 

so that the space-time wave function reads 
*l>(x,t) = (x\^(t)} 

N 

= £e-^*/ R C W)01*(0)>, (2) 

3=1 

with the the amplitude on site x denoted by the 
elementary trigonometric function — (x\j) = 

I 2 ■ jnx 
V N+l hm N+l ■ 
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Fig. 1. The space-time (x-t) representation of various initial 
wave packets moving ballistically in a chain. The classical time is 
T c i h(N + 1) and quantum memory effects are seen to develop 
with different initial wave packets giving different pictures, (a) For 
a local initial state ^(x, t = 0) = 8 x ,l in position 1 with N = 1000 
we see almost classical motion with recurrence at t = 1001. (b) For 
a spatially uniform initial superposition *$>(x,t = 0) = 1/VN in a 
chain with N = 100 we observe interesting interference effects, (c) 
For an initial uniform superposition of all stationary eigenstates 
with |^(0)) = ajb')j w ' tn a j = for any j, we have 

*(a;,t = 0) = y/2/N(N + 1) J^li sin amplitude at site x 

for a chain with N = 1000. The result resembles that of (a) where 
a wave is released at one end of the chain since in the spatial basis 
this initial state has higher amplitude near the end of the chain 
with diminishing oscillations, (d) For a Gaussian in space initial 
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superposition ^/(x,t = 0) 
patterns are also seen. 



; for N = 100 interesting interference 



The ballistic evolution is considered for various ini- 
tial superpositions ^(x, 0) = (a;|*(0)), chosen from 
expansions either in the spatial basis with |^(0)) = 
^lx=i a x\ x )i Ei=i \ a x\ 2 = 1 or in the eigenstate basis 
with |tf(0)> = Eti^lj), EtiKf = 1- The amph- 



tude Cx^ is real and the transformation from the spatial 
| a;) to the eigenstate \j) basis or its inverse from \j) to \x) 
are the same. In this case the computations of the space- 
time wave function %?(x,t) at various times are particu- 
larly simple by replacing the eigenvalues Ej, eigenvector 
amplitudes c? x in Eq. (2). For initial superposition in 
the spatial basis the normalized space-time probability 
becomes 



P(x,t) = \H>(x,t)f 



N 



N 



(3) 



and for superposition in the eigenbasis the normalized 
space-time probability is 



P(x,t) = \*(x,t)\ i 



N 



-iEjt/hjj) 



(4) 



Firstly the initial wave was released on a local state 
at one end of the chain with |*(0)) = |1) =4> *(x,0) = 
(x\^(0)) = S Xt i which from Eq. (3) at later times gives 
the probability 



P(x,t) = \H!(x,tf 



N 

i=i 



-iEjt/h CO* Jj) 



(5) 



in terms of the eigenvalues Ej, eigenvectors \j) = 

J2x=i c ^\ x )- This gives an almost classical evolution 
of a point particle seen in Fig. 1(a) where the particle 
at time T c \ = h(N + 1) returns to its initial position. 
The quantum revival resembles a classical particle mov- 
ing with constant velocity which reflects at the bound- 
aries at t = T c i and returns to its original position 1. 
Our second initial choice was the spatial superposition in 
space 1^(0)) = E^Li a x\x) with uniform a x — l/VN for 
all x, which gives from Eq. (3) the normalized space-time 
probability 



P(x,t) = \*(x,t)\ 2 = - 



N N 



£* 



x' = l 



(6) 



The corresponding space-time structure which arises 
from the initial condition of Eq. (6) is shown in Fig. 1(b) 
where we can observe interesting quantum interference 



features. The third choice we have made was not spatial 
but a uniform expansion on all stationary eigenstates \ j) 
choosing the initial superposition |^(0)) = J2jLi a j\j)> 
with aj = l/VN for all j. If transformed to the spatial 
basis this choice has a maximum at one end of the chain 
with decaying oscillations. At later times the normalized 
space-time probability is 



P(x,t) = \*(x,t)\ 2 = - 



N 

£< 



-iEjt/hJj) 



(7) 



Finally a Gaussian spatial initial superposition was cho- 
sen of the type 1^/(0)) = Y^=i a 'A x )i centered in the 
middle of the chain N/2. In this case after normaliza- 
tion the probability amplitude can be obtained from Eq. 
(3) with la, | 2 oc e -(--A72) 2 /(2- 2 ) 7 £^ =i \ ax \2 = 1; drawn 

from a Gaussian with a 2 = 1. 



The obtained space-time structures for the various ini- 
tial from the second to the fourth conditions are shown 
in Fig. 1. The white denotes high amplitude and the am- 
plitude is lower in the darker. The spectacular patterns 
shown are known as quantum carpets [2] and arise from 
scattering at the hard walls at the ends. The quantum 
intereference is more prominent for the uniform initial 
condition (Fig. 1(b)) where the picture resembles more 
closely the quantum carpets of Berry [1] . We can also see 
quantum recurrences more clearly in Fig. 1(a), (c) where 
the initial wavepacket has maximum amplitude close to 
the end site 1. The quantum carpets shown in Fig. 1. 
repeat indefinitely at time T c i, but less and less less ac- 
curately each time. 



B. Fractal dimensions for the probability density 



We observe the quantum evolution for a uniform ini- 
tial wave function along the chain as a function of space 
freezing the time, as a function of time at a single site or 
varying space positions equal to time. The correspond- 
ing pictures are shown in Figure 2 where interesting self- 
affine fractal oscillations are seen. For the spatial dis- 
tribution our simple computations for the fractal dimen- 
sion coincide with the value D x — 1.5 predicted by Berry 
This value coincides with the displacement curve of clas- 
sical random walks which gives Brownian noise. Quan- 
tum "randomness" described by the fractal dimension D x 
arises from the linear Schrodinger's equation, only due to 
quantum interference from scattering at the boundaries. 
It must be noted these curves are not self-similar but 
sclf-affine fractals. The computed time fractal dimension 
also gives the predicted value D t — 1.75. The D t is higher 
than D x as it can be seen form the more noisy oscillating 
of the probability density as a function of time t. We also 
display a space-time function along the diagonal x = t 
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where the self-affine fluctuations are much smaller with 
D x=t = 1.25. 
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t 
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Fig. 2. Fractal noise of an initially uniformly distributed in 
space wave packet for a chain of length N = 10000. (a) The abso- 
lute of the probability amplitude squared \P(x,t)\ 2 as function of 
space x at fixed time to = 100000. (b) A function of time t at the 
fixed site xq = 5000. (c) A function of both space and time with 
x = t. The obtained fractal dimensions are denoted in the log-log 
insets on the plot. 

The capacity (box) dimension of the self-affine y = 



y(x) curves is computed by counting boxes of linear size 
b.T along the x-axis axis and of size b. a along the y-axis. 
It must be noted that a similarity dimension cannot be 
defined since for a self-affine curve scaling is different 
along the two axes. We have used integer time units 
so the minimum box width along the x-axis is taken 
r = 1 and use a high-resolution grid for the amplitude 
along the y-axis. This is required in order to obtain the 
"local" non-integer fractal dimension D of the self-affine 
curves, while larger values of a finally give the well known 
"global" fractal dimension D = 1 [17]. The selected val- 
ues of a were of order of 1CP 10 for both the space D x 
and time D t dimensions. For the space-time fractal a is 
of the order of 7 x 10 -8 . The box dimension is defined 
by the scaling equation N(b,a) ~ b~ D , where N is the 
number of boxes needed to cover the curve [17]. A lin- 
ear fit of the scaling equation is displayed in the log-log 
plots of the corresponding insets of Fig. 2 which gives D 
from the slope. As pointed out above the obtained re- 
sults for D x ,D t ,D x=t are in agreement with the values 
for the fractal dimensions predicted by Berry [1] . For the 
other initial conditions considered we find similar fractal 
dimensions although with errors arising from the small 
sizes (see Fig. 3). 




Fig. 3. Time fractals for other initial conditions in a chain of 
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length N = 1000. (a) . A Gaussian spatial superposition at a fixed 
site xq = 500 where we compute Dt ~ 1.80. (b) A uniform initial 
superposition in the eigenbasis where we find Dt ~ 1.46. 



C. "Return to the origin" probability 



The "survival" or "return to the origin" probability is 
defined from the squared absolute amplitude overlap 



P (t) = |(*(0)|*(^ 



(8) 



which can be computed from Eq. (1). For an initial su- 
perposition in the spatial basis |^(0)) = J2x=i a 'A x ) onc 
obtains 



N N N 



j — l x—l x' — l 



(9) 



For an initial superposition in the eigenstate basis 
1^(0)) = ^2jLi a j\j) we have the simpler expression in- 
dependent of eigenvectors 



Po(t) 



N 



a,| e 

N N 



a -iEjt/h 



= EEM 2 M 2 ^ ( ^' )tA (io) 
j=i j'=i 
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Fig. 4. The "return to the origin" probability Po(t). (a) For 
an initially uniform superposition in space and a chain of length 
N = 100. The quantum revivals occur at an estimated period 
T r oc h(N + l) 2 »J 0.08(iV + l) 2 . (b) For an initially uniform su- 
perposition of stationary eigenstates in a chain of length N = 1000. 
In the inset a presentation of the probability density P(x, t). The 
quantum revivals which appear with period T c ; oc h(N + 1). 

The quantum revivals can be seen in Fig. 4 where Po(t) 
is shown versus t in a log- log plot. Po(t) is periodic in 
the early stages of evolution which becomes less clear 
later. The long-time overall asymptotic algebraic de- 



cay P (t) 



oc 



of ballistic systems can be seen pre- 



dicted from properties of Bessel functions. Fractional 
revivals with smaller strength are also seen. In the eigen- 
basis the expansion of the states j around the center 

_ QE— 

of the wave packet j Ej = Ej + (j — + ^r(j — 

7) 2 ^f + - = E-+2TTh[ 



J 

U-j) 



+ 



(3-D 2 



dj '2!' 

...1 allows to 



Q2-j • -■ - 3 • — -L T cl 1 2!T t 

find the time scales T cl _= h{N + l)/sm(jn/(N + 1)), 
T r = h(N + If I \tt cos(Jtt / '(N + 1))), etc. Of course, if 
|^(0)) consists of only one eigenstate with j = j the sys- 
tem shows no evolution. When the wave packets consists 
of many states j near the band center E w then the 
dominant time scale is T c i oc h(N + 1) with T r = oo. 
When it consists of states j near the edges Ew±2 then 
T c i = oo so the dominant time scale is T r oc h(N + l) 2 . 
In Fig. 4(a) for a uniform expansion in the spatial basis 
we observe period proportional to T r since it correspond 
to eigenstates close to E w ±2. In Fig. 4(b) the uniform 
distribution in the eigenbasis emphasises states from the 
band center E w which is the mean energy for this 
initial superposition, so that we observe a period propor- 
tional to T c i. 



III. DIFFUSIVE QUANTUM EVOLUTION 

We have also considered a chain with non-zero 
site potential taking the quasi-periodic value e x = 
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2cos(2irxM/N), x = 1,2, ...,7V for M, N successive Fi- 
bonacci numbers so that the rational M/N approaches 
the inverse golden mean ratio. In this critical case one 
obtains diffusive time evolution with the mean-square- 
displacement growing linearly with time (x 2 ) oc t which 
can be compared with the ballistic case (x 2 ) oc t 2 studied 
in Chapter II. It must be stressed again that although 
diffusive this is quantum evolution and not classical. 



A. Space-time structures 




Fig. 5. The space-time (x — t) representations for quantum 
diffusive evolution in a quasipcriodic chain of length N = 987. (a) 
For an initially uniform distributed wavepacket in space as a func- 
tion of space x and time t. (b) For an initially Gaussian wave 



packet in space. 

The obtained space-time structures for the diffusive 
case and two spatial initial conditions are shown in Fig. 
5. Again, the white colour denotes high amplitude and 
the darker lower amplitude. The figures show spectacular 
patterns due to quantum intereference from scattering at 
the hard walls at the ends of the chains but no quantum 
revivals. We show the space-time evolution for a uniform 
(Fig. 5(a)) and Gaussian spatial wave packet (Fig. 5(b)). 
The quantum interference gives very different results in 
this case since the obtained quantum carpets display self- 
similarity superimposed on the self-affinity found for bal- 
listic evolution. For the Gaussian in space initial state 
Fig. 5(b) we see slow quantum evolution with the system 
not reaching the boundary within T c /. 



B. Fractal dimensions for the probability density 




200 400 600 800 1000 



t 

Fig. 6. Fractal noise for diffusion in a chain of length N = 987. 
(a) An initially uniform wave packet at time t = 10 as a function 
of space, (b) As function of time at site xq = 500. 
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We obtain similar to ballistic case self-affine fractal 
dimensions displayed in Fig. (6). Apart from self- 
affine the curves are also self-similar fractals with di- 
mensions defined in [16]. From the normalized probabil- 
ity measure P(x,t) with ^2 x= iP{x,t) = 1 we compute 

" E»=i P ( x ' *) ^P(x, t), lnJ2x=i P i x , and from the 
linear fits vs. In TV obtain the space fractal dimensions 
Df, D%, respectively. Our results gave give the entropy 
fractal dimension Df w 0.85 for the uniform and « 0.43 
for the Gauss while for the correlation dimension we ob- 
tain ss 0.77 for the uniform and s» 0.37 for the Gauss. 
It must be pointed out that the fractal dimensions D± , D2 
depend strongly on initial conditions (see Fig. 7). 




2 2.5 3 3.5 4 

log N 



Fig. 7. Entropy vs. logN for the diffusive model for an ini- 
tially (a) uniform and (b) gauss distribution, both in space. The 
information dimension Df s=s 0.85 for the uniform and Df w 0.43 
for the Gauss is obtained form the slopes of the linear fits. Similar 
results exist for the correlation dimension D2. 



C. "Return to the origin" probability 



The "return to the origin" probability now shows many 
fluctuations and if we define the time integrated correla- 
tion function C{t) = (l/t)J*P (t')dt' vs. t. We also ob- 
serve the long-time asymptotically algebraic overall de- 
cay Po(t) oc -4= known for diffusive systems. 




1 10 100 1000 10000 



t 

Fig. 8 . The "return to the origin probability" for an initially 
uniform distributed wavepacket in a chain of length N = 987. The 
inset shows the log-log scale of the correlation function C(t) vs. t 
with diffusive correlation behavior. 



IV. DISCUSSION-CONCLUSIONS 

We studied the milder when compared to classical 
quantum ballistic and diffusive evolution of wave pack- 
ets. It is different from that of classical point particles, 
nevertheless, shows fractal noise in the form of self-affine 
curves. Our calculations are done by discretizing the 
space into a tight-binding chain lattice. For the bal- 
listic case semiclassical diffusive evolution does not ex- 
ist for small t, to cross over to ballictic for longer times 
[18], since the ballistic law is satisfied for all t. In this 
case we display quantum memory effects by observing 
different evolutions for various initial superposition wave 
packets ^(x, t). However, in agreement with Berry's find- 
ings for a "particle in box" a self-affine fractal space-time 
P(x,t) — \^(x,t)\ 2 is found. For diffusive systems the 
evolution apart from self-affine is also self-similar with 
no recurrences. 

The quantum evolution could be efficiently imple- 
mented by a quantum computer so that one can use 
the simulation to solve certain computational tasks by 
creating more efficient algorithms. In fact, Po{t) can re- 
veal the prime factors of an encoded number [12]. On 
the other hand, fractals known from non-linear physics 
also appear in quantum physics, despite the fact that 
quantum mechanics is linear. The space fractal dimen- 
sion D x = 1.5 obtained for the spatial distribution cor- 
responds to Brownian l// 2 noise familiar from random 
walks. Therefore, techniques can be transferred from 
non-linear physics to quantum mechanics and vice versa. 
Finally, quantum intereference can make itself obvious 
in scattering from nanostructures, such as carbon nan- 
otubes. It can possibly influence the temperature inde- 
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pendent white noise known as "shot noise" for measure- 
ments of current in nanostructures [19]. 

In summary, classical systems are diffusive due to ei- 
ther non-linear effects, which may lead to deterministic 
chaos, or introduced randomness from outside (flipping 
a coin) as for the random walk. In these cases irre- 
versibility naturally arises with convergence to a Gaus- 
sian steady state distribution for any initial state. In 
quantum systems the evolution is smoother, wave-like, 
reversible, with no convergence to a steady state, since 
quantum memory develops. However, we show for simple 
lattice systems fractal noise due to quantum interefer- 
ence from scattering at the boundaries. The noisy quan- 
tum evolution seems not only independent of chaos but 
is present for all kinds of evolution. 
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